TECHNICAL MEMORANDUMS 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



No. 829 



METHOD OF CURVED MODELS AND ITS APPLICATION TO 
THE STUDY OP CURVILINEAR FLIGHT OF AIRSHIPS 

PART I 
By G-. A. G-our jienko 

Central Aero-Hydrodynamical Institute, Moscow, 1934 

^ ^ 



Washington 
June 1937 



NASA TKftnlcal Ubrvy 




N.A.C.A. Technical Memorandum Nov 829 
PREFACE 

The idea of the desirability of testing curved models 
for studying the 'aerodynamics of airplanes in curvilinear 
flight (for example, in a loop), I have expressed in my 
lectures on the dynamics of flight a very long time ago r 
as early as 1918, or even earlier, hut I did not intend 
to publish it.* And this idea did not materialize until 
Comrade Gourjienko devoted ahout a year to the treatment of 
this question. 

The possibility of substituting for one phenomenon - 
the flight of a straight model along a circular path - 
another phenomenon, the flight of a curved model along a 
rectilinear path, is based on the following proposition: 

If two flows under the influence of bodies immersed 
in them experience the same geometrical distortion of the 
lines of flow, the forces acting on the bodies may be re- 
garded as approximately equal. 

In simple cases, the correctness of this assumption 
is full?/ - apparent. For example, the rectilinear motion 
of a flat plate and the circular motion of a plate bent 
along the circumferential arc, along which the motion oc- 
curs (figs. 1 and 2). 

The analogy holds true if in place of a thin plate, 
we use a symmetrical streamlined form (the so-called 
"rudder- form" ) , as is shown in figures 3 and 4. 

Passing to theoretical forms of wings, consisting of 
a basic arc, bent, along a circumferential. .arc* and enclosed 
by a streamlined form, for example, the inversion of a 
parabola or an ellipse, we see that the theoretical lifting 
force of such wings in plane-parallel flow is algebraically 
composed of two components: the first depending on the an- 
gle of attack, and the second depending on the concavity 
of the basic arc and proportional to it. 

For the small arcs of one span, which differ according 
to their concavity, the second component of the lifting 
force is proportional to the curvature of the basic arc. 

♦Since I occupied myself with aerodynamic questions only 
in conjunction: with other questions (reliability of design 
[designed stability ?], dynamics), and the ideas that oc- 
curred to me I expressed in lectures or discussions, in 
order that persons, working in this particular field, 
might make use of them. 
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Causing the curved theoretical wing to move under a 
zero angle of attack in three ways: (l) rect il inear ly , 
(2) along the circumferential arc, coinciding with the 
"basic arc, and (3) along the circumferential arc of the 
somewhat larger radius, we find in the first case (fig. 
5) a lifting force proportional to the curvature of the 
basic arc; in the second case (fig. 4), we find a lifting 
force equal to zero, and in the third case, the intermedi- 
ate case (fig. 6), we can expect no other lifting force '•• 
than the one which is proportional to the difference be- 
tween the curved "basic arc and the circumferential arc 
along which the wing moves. 

On figure 6, 0 is the center of curvature of the 
basic arc; 0 X is the center around which the wing moves. 

The continuous arc AB is a basic arc, with its cen- 
ter at 0« 

The dotted arc AB is the path of the points A and 
B when they travel around the center 0; the straight 
line (dotted line) is the path of the same points when 
they travel along the straight line. 

Finally, if the wing moves under an angle of attack, 
we may here too extend the theorem concerning the lifting 
force to the case of motion along the circumferential arc, 
the angle of attack denoting the angle of turn of the wing 
with respect to its position. AB (figs. 5 and 6). 

With the considerations cited, the basic assumption 
for the method of curved models can be regarded, if not 
as demonstrated, then as entirely supported by evidence, 
at any rate for the cases when the radius of curvature of 
the wing trajectory is comparatively large with reference 
to the length of the chord of the wing. 

Of course, the curvature of the path involves correc- 
tions which must be regarded more appreciable when the 
lifting force of the wing is directed away from the center 
of rotation, and less perceptible when the lift force is 
directed toward the center of the trajectory curvature. 
In aviation and in aerostatics, exactly the latter case is 
involved, which is more favorable for the application of 
the elementary theory. 

It seems to me that for plane-parallel flow we can 
also give an exact hydrodynamic theory for the wing in 
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circular motion, using as a "basis the writings of P. A. 
Walter, on the theory of rotating turbine "blades in a 
converging flow. 

As regards the motion of three-dimensional "bodies, 
analyzed "by the author (dirigibles , airplane fuselages, 
air "bombs), we must not expect here to obtain a hydrody- 
namic streamline theory, since such a theory is absent 
even in the simplest cases of rectilinear motion, of a 
straight dirigible under an angle of attack. Our analogy 
is fully applicable also to this case of motion. 

And, we may expect that the experiments and rechecks, 
according to the method of curved models, will outdistance 
the theory. • 

It should be noted here that the principle of the ap- 
parent curvature of the wing was applied by Glauert when 
he calculated the effect of the rotation of the wing on 
its lift. Having much less theoretical evidence, he never 
theless arrives at results that are in good agreement with 
experimental results. 

In substance, the Glauert method is as follows: If a 
flat plate (basic arc of a symmetrical wing) moves for- 
ward and rotates upon going into a dive, the geometrical 
angles of attack (with relation to the flow unaltered by 
the wing) will be different for the different elements of 
the basic arc (fig. 7). 

In order to produce. such geometrical angles of attack 
for the wing in translatory motion, we must bend its basic 
arc as shown in figure 8. 

On figure 7 the arrows A 'A, b'b, c'c, ..... k»k, 
B'B show the direction of the lines of the unaltered flow 
with reference to the various elements of the plate which 
is moving to the left and is rotating. 

Figure 8 shows a curved arc, which the lines of the 
unaltered flow approach at the same angles as in figure 7. 

This same analogy can be established in another way» 
Let the plate, which is curved along the arc of the circle 
move in still air along its direction, as in figure 2, at 
a linear velocity u and at an angular velocity: 
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The angle of attack of each element will he equal to 

zero. 

Now let us stop the rotation of the plate and let us 
cause it to move along the tangent to the circumference. 
The plate will then assume the position shown in figure 8. 
But this will last one moment, and while it moves along 
the circumference, when the rotation of the plate ceases: 
this is analogous to the application to it of the angular 
velocity a) in. the reverse direction. But, while main- 
taining the movement of the center of gravity, we cannot 
ohtain the conditions of figure 8 for any length of time. 

However, the application of the theory of the "dynam- 
ic camber" to the oscillating wing gives very good results, 
as is found "by tests on vibrations of wings in an air flow, 
"by calculation of the damping effect of the wings and "by 
determination of the rotary derivatives "by means of the os- 
cillation method. 

We should expect the method of curved models to "be 
much more exact, from the theoretical point of view, not 
to mention the simplicity of testing a fixed model in the 
tunnel. 

And we should "be grateful to Comrade G-ourjienko for 
his initiative and lahor in working out this method and 
making it practicable. 

Professor Y. P. Yetchinkin. 
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METHOD OE CURVED MODELS AND ITS APPLICATION TO 
THE STUDY OP CURVILINEAR FLIGHT OP AIRSHIPS* 
By Cr. A. G-ourjienko 



INTRODUCTION 



Up to the present time, the experimental study of the 
aerodynamics of curvilinear flight and of the very closely 
related question concerning the dynamic stability of air- 
ships as well as other aircraft has "been in the embryonic 
state. The reason therefor lies in the fact- that in solv- 
ing each of the questions mentioned, in addition to the 
usually given data on aerodynamic resistances (forces and 
moments) in rectilinear flight, that can he sufficiently 
determined by the generally recognized method of testing 
in wind tunnels, it is also necessary to include the forces 
and moments due to the presence of angular velocity. 

In order to determine these forces and moments, usu- 
ally expressed by the so-called resistance rotary deriva- 
tives, special equipment was hitherto required. 

One of the most accurate (in the kinematic and dynam- 
ic sense) methods was the use of whirling arms, where there 
was reproduced to some extent the phenomenon of circular 
flight. However, in practice this method is not very suit- 
able for wide application, on account of the following 
reasons: (a) the extreme difficulty in obtaining exact 
measurement of the aerodynamical forces and moments acting 
on the model, taking into account the inevitablo influence 
of tho centrifugal loads; (b) the complicated air stream 
that is created during the motion; (c) tho rather compli- 
cated and expensive equipment. 

Another generally recognized method is the small os- 
cillation method. This method does not require any large 
special installation, and allows experiments to be carried 
on in ordinary wind tunnels - although a special, occa- 
sionally rather complicated apparatus, is used. 



♦Report No. 182, of the Central Aero-Hydrodynamical In- 
stitute, Moscow, 1934. 
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We shall not dwell upon the well-known principle un- 
derlying this method, hut will merely point out that by 
means of extremely laborious experiments and working up 
of data, it enables us to determine n nlv the rotary deriv- 
atives of the moment. Besides, the accuracy of the deter- 
mination of this derivative is very low, since in addition 
to the large number of different - only approximately cor- 
rect - propositions, lying at the very basis of the theory 
of the method, this derivative, in practice, has to be ob- 
tained by measuring an extremely large number of quanti- 
ties, the majority of which has no connection with aerody- 
namics, as, for example, the moment of inertia, friction 
in the apparatus, rigidity of the springs in the equips 
ment , etc. 

The inevitable errors in the measurement of all these 
"ballast" quantities, as well as the damping effect itself, 
produce very large, errors in the final result. 

As we have already indicated, we can obtain with this 
method only the rotary derivative of the moment. The ro- 
tary derivative of the lateral (or lifting) force can be 
determined from the derivative of the moment only by means 
of various approximate repeated calculations (English 
method) . 

But the rotary derivative of the drag, which is ex- 
tremely important, for instance, for investigating the de- 
celeration while turning, cannot be determined at all by 
the aid of the method of small oscillations. 

Moreover, neither of the existing methods in any way 
answers the very important question concerning the distri- 
bution of aerodynamical loads on the various parts of air- 
craft (hull and tail surfaces of a dirigible, wings of an 
airplane) while turning and during curvilinear flight. 

The practical impossibility of experimenting on the 
distribution of pressure on whirling arms is due to the very 
complicated- arrangement of the pressure tubing on the bat- 
tery gage with the revolving model, and is due to the com- 
plicated calculation of the centrifugal forces, acting on 
the columns of air in the tubing. 

In view of the above, it was decided in the aeronauti- 
cal section of the Experimental Aerodynamical Department of 
the Central Aero-Hydrodynamical Institute, upon the initia- 
tive of Professor K. K. Fed;iaevsky , to work out as complete- 
ly as possible both theoretically and experimentally an en- 
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fly— arg—p-&-3-a -ib-l-e— b-o^h — t^^o.r^.tA"e&l^^~a«'d^exp.e.r^i-m^tt-frarl"l'y te sn 
entirely new aerodynamic principle of testing curved mod- 
els. Several allusions to the possibility of applying 
this principle were made in foreign literature dealing 
with dirigible aerodynamics. As far as. we know, Professor 
V. P. Vetchinkin (see Preface) was the first to speak of 
the possibility of applying this method. Put, unfortu- 
nately, up to the present, no one (either abroad or here) 
has developed this method. Therefore, we had to work out 
this method entirely independently'. The work w a s done in 
191*3-34, and was financed" by NIO of the dirigible factory 
(Scientific Research Division of the dirigible factory). 

The basic premise of the method is as follows: In 
the ordinary wind tunnel with rectilinear flow it is nec- 
essary, for the model installed in it, to create a kine- * 
matically similar pattern of flight along a circular tra- $ 
jectory. This can be obtained, in the first place, by \0^^e" 6 ' 
bending the model in a special way and, secondly, by cre*A 
ating across the tunnel a constant velocity gradient J 
(since, during flight along a circle the velocity of mo- 
tion of each point of the dirigible is a linear function 
of the radius of its trajectory). Thus, by testing the c , 
curved model on the conventional aerodynamic balance, we ^ 
obtain the forces and moments acting on the model in the p J 1 " 

proposed motion along a circle. This will enable us to r & 
obtain the rotary derivatives of the moment, of the lat- f A^,} .f 
eral force and of the drag completely independently, withy ^ ' 
accuracy and simplicity, by means of typical, simple ex- v*j> 1 
periments in the tunnel. Moreover, the simple arrangement 
of the pressure tubes of the curved model will enable us 
to study the distribution of pressure in curvilinear 
flight. 

In the first part of this paper we shall, present the 
theoretical side of the problem of constructing curved mod- 
el forms and the method of testing the model. In the sec- 
ond part we shall present a detailed account of the first 
experiments according to the given method, carried out with 
a curved model of the nonrigid airship V-2, and a compari- 
son of the experimental results with some data of full- 
scale flight tests made with this airship in 1933, 

The author expresses his sincere appreciation to Pro- 
fessor K, K, Pediaevsky for his valuable and highly compe- 
tent advice and assistance in the solution of various prob- 
lems. 
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PART I 

1. POSSIBILITIES OP ROTATING- CIRCULAR MOTION 



Before we come to the theory of constructing curved 
model forms, i.e., to the solution, as we shall see, of a 
purely geometrical problem, it is necessary to discuss the 
possibility of rotating circular motion, in a manner simi- 
lar to that applied in rotating rectilinear motion in ex- 
perimental aerodynamics. 

We must admit that from the point of view of theoret-- 
ical hydrodynamics such a rotation of motion is inadmissi- 
ble. While in moving along a circle the dirigible encoun- 
ters a mass of absolutely calm air, when the motion rotates 
(if we make the air rotate like a solid body)* the flow be- 
comes vortical at every point in space, the intensity of 
the vortex being X = 2iUy , where ujy is the angular ve- 
locity of rotation. This can be easily verified by tak- 
ing the value of the velocity circulation along any closed 
contour within the indicated hypothetical flow. If we 
were to consider the given assumption as an obstacle when 
passing from circular flight to a stationary model in cir- 
cular flow and thence also to the curved model, it would 
be necessary to question also the validity of the use of 
the wind tunnel altogether - in which, as is well known, 
the flow is vortical, nonpotential , and nonrectilinear . 

We may assume that the errors, due to the fact that 
rotation of circular motion does not sufficiently comply 
with natural law, are not larger than those which are ob- 
tained by means of the usual experiments in the tunnels. 

Consequently, the method of curved models, although 
it is not absolutely exact, is nevertheless a close approx- 
imation to reality. 

2. CONSTRUCTION OP THE CURVED MODEL PROFILE 

As an initial principle for the construction of the 
curved model profile, let us make the following indis-^ 



*Incidentally , it should be mentioned that such motion of 
a fluid is physically impossible; that is to say, such a 
flow can exist for only an infinitesimal given space of 
time. 
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putable assumption: If every, small section of the sur- s$ - 

face of the curved model forms - with the direction of Lfty ^fi 

the" rectilinear 'f low in the wind tunnel the same inter- itff 



secting- angle as will be formed hy a corresponding sec- *f' 
tion of a straight model with the direction of its linear 
speed while turning along a circle, the kinematic, and 
consequently, the dynamic similarity will he observed. 

M 

The similarity in the sense of the magnitude of. ve- .pJV . 
locity of motion of the sections considered is maintained H : \ff 
hy creating in the tunnel the constant velocity gradient 
referred to above. However, it does not seem possible to 
solve the problem in such complete form. It is not diffi^ 
cult to observe that by curving the model, the length of 
the arc of the meridional contour inevitably becomes some- 
what shorter on the windward side and somewhat longer on 
the lee side. To introduce this degree of contraction and 
elongation into the solution of the problem is impossible, 
because, in order to determine this degree, it is necessary 
to have the problem already solved. 

However, the problem can be solved with a very great 
degree of approximation if, instead of the surface, we ex- 
amine the kinematic similarity of the motion of the axis 
of the model and of an infinitely large number of cross 
sections. After making all the calculations concerning 
the axis and the cross sections, we shall try to find the 
numerical value of that error, which results from substi- 
tuting for a complete solution an approximate solution, 
and which, as we shall see, also is due to the indicated 
contraction and elongation of the arc of the meridional 
contour. 



And so, we 'shall proceed to derive the equations of 
the curved axis and the curved cross sections. 

Let us assume a dirigible turning along a circle with 
radius R 0 (fig. 9). Considering that its center of vol- 
ume 0 moves along the circular trajectory A - A with a 
velocity u 0 , the theory and the experiment show that the 
axis of the airship will form an angle of attack 3 Q with 
the direction of velocity v 0 . 

Besides, the rudder is deviated, as shown on the dia- 
gram, at an angle 8. In addition to the basic trajectory 
A - A, the diagram also shows a series of other trajec- 
tories along which the points of the airship axis move. 
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From the center of rotation 0 let us drop on the 
axis of the model the perpendicular OB. Then the angle 
COB = 8 0 , as an angle with the sides perpendicular to 
U 0 C and BO. From the triangle BOC, we obtain: 

OB = R 0 cos P 0 and BC = x 1q = R 0 sin B 0 

How let us take on the airship axis the arbitrary- 
point D, which moves during rotation along the trajec- 
tory m - m; that is, concentric with relation to A - A. 
Let us connect point D with 0 and draw at D the tan- 
gent to m - m. This tangent gives the direction of veloc- 
ity u of point D. The angle uDB will be the local 
angle of attack S of the infinitely small segment of the 
airship axis at point D. 

From the diagram we see that 8 > S 0 . 

If along the axis we take point D closer to point 
B than is shown on the diagram, angle 8 begins to get 
smaller and becomes 0 at point B. 

Designating the variable distance BD by x x , we 
find 8 = f(z 1 ). From the triangle BDO , where angle 
BOD = 8, we obtain: 

x x = R 0 cos 8 0 tan 8 

whence 

tan 8 = 



*1 



R 0 cos 8 0 

and 8 = arc tan (l) 

R 0 cos Po 

it is clear that when x^^ = x x , 8 = 3 0 . 

How, for the circular motion, let us substitute recti* 
linear motion, shown in figure 10. 

Starting from the basic principle, it is necessary to 
bend the axis of the model in the direction indicated on 
the diagram, in such a way that the angle of attack along 
the axis of the model changes according to the same law 
(l) as in the case of curvilinear flow. 
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Naturally, when the axis and the cross sections are 
curved, it is necessary to preserve their lengths. This 
seems all the more necessary since, "by preserving these, 
lengths, we maintain the distribution and magnitude of the. 
velocity circulation along the axis and along each of the 
cross sections, which is very essential for the similarity. 

Thus, x x will denote the length of the arc of the 
curved axis from B to D. 

Now, taking on figure 10 point B as the origin of the 
coordinates x, z',* the direction of which is indicated, 
we obtain the equation of the curved axis. 

The angle of inclination of the tangent to the curve 
sought at point D must he equal to (3, whence 

tan 0 = (2) 
dx 

On the other hand, it has already been shown that 



tan 0 = 



R 0 cos P 0 



where x x is the length of the arc, and is equal to 

x 



0 

Then the expression (2) is transformed as follows: 

x 

&*± = i r n ~ TlzSf dx (s) 

dx S 0 cos 3 0 / v \dx / 

o 

The equation obtained is also the equation of the 
curved axis of the model. Let us solve it. For this pur- 
pose, let us differentiate both sides and, for the sake of 

brevity, let us take — 1 = n: 

E 0 cos 8 0 



dx 3 <J \dx / 



♦The direction z 1 coincides with the direction z, taken 
as a reference. Thus, the origin of the coordinates is at 
a distance of z' c from the axis of the tunnel, that is 
to say , z = z • - z 1 c . 
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dz ' 

let us substitute the new variable u = — — . Then 

ax 



& = n A + u 2 
dx v 



Dividing the variables, we obtain: 

du 



dx 



n Jl + u £ 
Let us integrate: 



II" du 
x = — 



The integral obtained is the "long logarithm" and is 
expressed: 



=« ^ In (u + /l + u a ) + 0 (4) 
n v 



The constant 0 we determine from the initial condi- 
tions. When x-i = 0, 

x = 0 and tan 8 = = u = 0 

dx 

Substituting, we find that C = 0. 

Let us solve the equation obtained with reference to 
u, for which purpose we find from (4) 



and 



u + Jl + u 2 = e n3C 



nx 
e - u 



Jl + u 3 
Squaring both sides, we get: 

1 + u a = e 2ns - Sue 11 * + u 3 
We reduce by u s and find u = f(x): 



(5) 
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Substituting u = , we get: 

dzj _ e a nx - 1 
dx 2e nx 

dividing the variables and integrating, ire find: 



r e anx - 1 



2e 



nx 



dx 



(6) 



Performing division under the integral and dividing 
the right side into two integrals, we get: 



e nx dx - i 



e~ nx dx 



Performing integration and substituting, we get; 



.! = -"L [e nx - 1 + e" nx - 1] 



e nx + e -nx 



- 1 ! 



e nx + q— nx 

In analysis, the expression tj , is called a 

hyperbolical cosine and is denoted by ch nx. Thus 

z« = i [ch nx - 1] 
n 



Taking n = 



R 0 cos P Q 



we find that 



= R 0 cos Po (ch ir -f--p- - l) 



3c 

This is also the equation of the curved model axis, 
ing up the parenthesis, we obtain: 



(?) 
Open- 



z» = E, 



cos p 0 ch ( 



x 



R 0 cos P 0 / 



R Q cos P 0 



This is no other than a "catenary" with the parameter 
R 0 cos P Q = OB (fig. 9) with the origin of the coordi- 

#S ***** Mr* 4.***),*** ****** 

<*4»*j cfr-c ar^r &4*sS M cwrS ****** no * 
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nates transferred to the point of the minimum 
form is assumed by a perfectly flexible thread 
at two points and subjected to its own weight, 
proportional to the length of the thread. The 
tained is very convenient for constructing the 
for hyperbolical functions there are available 
bles. (See Hutte, vol. I; or Keiichi Eayashi , 
mehrstellige Taf eln der Kreis- und Hyperbelfunktionen , " 
[Tables of seven or more places of circle functions and 
perbolical functions], Berlin, 1926.) 



z 1 i Such a 

suspended 
which is 
formula ob- 
curve, since 
detailed ta- 
"Sieben- 



und 



hy- 



Generally , the ordinates of the meridional contour of 
the dirigible and other data on the hull are given as func- 
tions of the distance along the axis of the dirigible from 
its nose. Therefore, it is useful to have the expression 
which enables us to find the quantity x for the construc- 
tion of the curve according to formula (7), given the quan- 
tity xi, that can be measured off along the arc of the cat- 
enary from the origin B of the coordinates. 

The relation x = f(x x ) is readily obtained by inte- 
grating in the expression (2a): 

x 



x, = 




1 + 



dx 



Above there was substituted 
ing to formula (5), was equal to: 



d zi _ 
dx 



u, which, accord- 



u = 



3 snx _ 1 e nx _ e -nx 



2e 



nx 



This expression is no other than the hyperbolical 
sine of the argument nx; that is, to say, 

u = sh nx 



Substituting in formula (2a), we get: 

r-. X 



*1 



1 + sh nx dx 



According to the basic relation between the hyperbol- 
ical functions of one argument : 



ch 3 nx - sh a nx = 1 



N.A.C.A. Technical 



Memorandum No. 829 



11 



we 



find that 



1 + sh nx 



sh nx 



= 1, 



wherefore , 



x 




ch nx dx = — 



n 



L 



sh nx 




— sh nx 



(when x = 0 sh nx = 0 ) . 

Solving the expression obtained with reference to x 
we have the relation sought: 



(The symbol Ar denotes the inverse hyperbolical function 
and comes from the word "Area.") 

Now let us proceed to the study of the curvature of 
the cross sections. The method of investigation is the ■ 
same: We shall examine the variation of the angle of in- 
tersection along the horizontal diameter of each cross sec- 
tion. 

The direction of one of these diameters, passing 
through the arbitrary point D' on the axis, is represented 
on figure 9 by the straight line a - a. 

let us derive the equation of the curve a - a (fig. 
10), which changes from the straight line a - a in pass- 
ing from curvilinear flow to rectilinear flow. 

It is obvious that the curve sought will also be a 
catenary, normal at point D to the curved axis, since the 
variation of the angies of attack along a - a on figure 
9 does not differ theoretically in any way from the case 
examined in connection with the axis of the model. Let us 
prove it. 

Taking on the straight line a - a any arbitrary 
point Dj, , let us construct for it the local angle of at- 
tack designated by $ . Let us connect point Dj, with 0 
and drop from 0 the perpendicular 0'B l on the direction 
a— a. Then the angle' D 1 0B X , which is equal to $ , is 
expressed as follows: 




(8) 
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OBj. Xi 

whence i = arc tan — — - - (9) 

Xl 

The expression ohtained is completely analogous to 
the formula (1) derived for the axis of the model. The 
distance D 1 Bj coordinates the distance of point Di 
from Bi , similarly to the manner in which in formula (1) 
the distance coordinated point D with reference to 

B. The quantity x x (formula (9)) for the entire straight 
line a - a plays the same role as the quantity R 0 cos p 0 
in formula (1 ) •■ 

It is obvious that, for the curve sought, point B x 
will perform functions analogous to the functions per- 
formed by point B for the curved axis. At Bi the an- 
gle ■b will be equal to zero; that is to say, B x will 
be the origin of the coordinates of the catenary sought. 

In order to derive the equation for the curve sought, 
let us turn to figure 11, It shows the curved axis BB 
in terms of the coordinates x, z', and the curve a - a 
sought. From what we have said in making an analogy with 
the above, it follows that, for the curve sought, we should 
'take B x as the origin of the coordinates, and draw the 
axes of the coordinates parallel to the axes x, z'. De- 
noting these axes by r and t, let us write, analogously 
to formula (7), the equation for the catenary a - a in 
terms of the coordinates r, t: 



t = x x (ch - r - - l) (10) 



Xi 

It is obvious that for the line a - a, xj.= constant 
When passing from one cross section to another, that is to 
say, when varying x x , we naturally change the position 
of the origin Bj. of the coordinates. As x x decreases, 
point Bj will travel along the curve BiBj" , shown on 
figure 11 as a dotted line. Moreover, the catenary will 
"contract," assuming for a certain point B x ' the shape 
a' - a', indicated on the diagram. This_J.s due to the 
fact that the arc of the catenary (10) B X D, which is 
equal in length to the distance 

B7D = H 0 cos Po = n 
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(fig. 9), does not change when x x changes; that is to 
say, . 

B X D = Bi'D' = Bj'B = R 0 cos 6 0 = J = constant (11) 

When Xi = 0, the catenary will contract into a 
straight line, disposed along the axis z'. On crossing 
the axis z 1 , the catenary will "open up" again. 

Naturally, all the chords of the cross sections, par- 
allel to the horizontal diameter, will deflect along the 
segments of the same catenaries as the horizontal diame- 
ters, and all the vertical chords and diameters will re- 
main straight. This is due to the fact that the motions 
examined are plane-parallel. Thus, the cross sections will 
curve along the cylindrical surfaces. 

Now, let us express the whole family of catenaries in 
terms of the coordinates x, z'; this is very convenient 
in constructing the profile of the meridional contour of 
the curved model. 

Let us denote the coordinates of point D "by xj> 
and z'j), and the coordinates of point B 1 "by x 0 , z 0 1 . 
Then the formulas for the change of coordinates will he 
written in the following form: 




(12) 



Let us find the values x 0 and z 0 1 . Let us desig- 
nate the length of the arc B,D hy r, . According to 

D 

condition (11), we find that 

ri D = n 

On the other hand, the projection of the arc Bj.D on 
the axis r, equal to r-p , analogously to equation (3) 
for the catenary a - a, is written: 

ri B 1 
. r D = x i Ar sh -— = x x Ar sh n5E^ 

We may disregard the symbols r^ and r x , since we 



14 



N.A.C.A. Technical Memorandum No. 829 



are here interested only in the absolute values of the 
segments and arcs. 

From figure 11, we see that 



whence 



x 0 = x D + x x Ar sh 



(13) 



In order to obtain z 0 ' , let us find the ordinate 
t 0 of point D. Considering that 

r-Q = x, Ar sh — — 



nx. 



according to formula (10), we find: 

t B = x x (ch |2 - l) = Xl ch (Ar sh ^i-) - 1 



t i on s 



According to the basic formula for hyperbolic func- 



ch 8 cp - sh 2 cp = 1 



the preceding expression may be written as follows: 



1 + sh" Ar sh 



or 




/l. 



= n </( nx i 



) + 1 - X! 



But, starting from the equation of the curved axis 
z'j, = \ [ch nx D - 1] 

in which, according to expression (8), we may substitute 

nxji = Ar sh nx x 

we find that 

Z, D = n ch ( Ar sh nx i> ~ i = n / (nx!) 3 + 1 - J 
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Thus, t 



D 



= J /(n Xl ) S + 1 - Xl = z » D + J 



Prom figure 11, we see that z 0 1 = z - t D . Conse- 
quently, 

z Q t = z' D - z« D - J + xj = x x - I (14) 

How, knowing x 0 and z 0 » as functions of the quan- 
tities xi and x D , that is to say, of the quantities 
coordinating point D in the x, z 1 system, let us repre- 
sent the formula for the change of the coordinates (12) -as . 
follows : 



r = x - xj, - Xi Ar sh- — 
t = z* - x, + A 

* n 



1_ 
1 



(15) 



Then the catenary a - a (formula (10)) in terms of 
the coordinates x,z' is expressed in the form: 



Z 1 - X n + — = Xi 

1 n A 



Finding from here z r = f(x) and taking x x out of the 
parenthesis, we obtain: 

z« = Xi \ch(-——2 - Ar sh -J (16) 

This is the equation of the family of catenaries 
sought. 

It was shown above that each catenary of the family 
(16) must intersect the line of the curved axis at right 
angles. Let us verify this circumstance. 

Prom analytical geometry we know that for mutual or- 
thogonality of two curves, it is necessary and sufficient 
for the values of the derivatives of these curves at the • 
intersecting point to he inverse in magnitude and sign. 
Let us find the derivative with respect to x from formula 
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(16):* ■ 

|zi = . Xl sh . Ar sh JJ) JL = sh (f^B - Ar S h 

dx \ x x nx x y Xj \ x x nij^^ 

The -value of the derivative at point D is obtained 
"by substituting x = x^; we have: 

d- 1 ) - * h C~ Ap sh = - -i- (18) 

(x=x) D 

According to the preceding, the value of the deriva- 
tive from the equation of the curved axis will be: 

l-S- = u = sh nx 
ox 

substituting for x its value, expressed as Xx accord- 
ing to formula (8), we find that 

|zl = sh [ Ar sh nXi ] = nac (19) 
ox 

A comparison of the results (18) and (19) shows that 
the catenaries (7) and (16) are orthogonal to each other 
at point D. 

Now, it remains for us to find, analogously to the 
above, the expression which will enable us to lay off on 
the curves a - a from point D the magnitude of the lo- 
cal radius of the cross section of the model. Denoting 
this local radius by y x , let us assume (fig. 11) that it 

is plotted along the curve a - a from point D above 
and below; that is to say, 

DA = DB = y x 

Thus, points A and B belong to the meridional con- 
tour of the model. Since it is difficult to plot along 
the arc, let us try to find the abscissas of the points 
A and B in terms of the coordinates x, z' for a given 
y 1 . Designating these abscissas by a and b, we find 
the length of the arc AB of the curve a - a: 



*Prom the theory of hyperbolic functions we know that 
ch cp = sh cp and sh (-cp) = -shcp 

acp 
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.2 

dx 



AB = 2y x - [ J 1 + (tr) 

a a i 

Introducing the value of the derivative -~- accord- 
ing to formula (17) and substituting according to formula 

J 1 + sh 8 cp = ch qp . 

we obtain: 

. — . n ✓x - x D i \ 

AB = 2 7l =y ch ~- - Ar sh — ) dx 

Integrating and considering that for the given curve 
a - a, xj) = constant and x x = constant, we get: 

AB = 2 n = Xl [ sh (—2 - Ar sh - 

- sh (i-TJS - Ar sh £)] (20) 

How let b = xp, that is to say, let us assume that 
therejvas integrated not the segment of the arcs AB , 
but AD = y l . Then, substituting b = x D in the preced- 
ing equation, we find: 

15 = Vl = - x Tsh - Ar sh ~L) + -L.] 

{_ V xi nxi/ nxj.J 

Prom the equation obtained, we determine the abscissa a 
sought : 

a = x D + x, [ Ar sh -~ - Ar sh (|± + —)] (21) 

Letting, in equation (20), a = x D , that is to say, 
assuming that arc DB = y x was integrated: 



D3 
whence 



= 7x = *i [sH (--57— - ^ s h Hx~) + Hx^ J 



b = x D + x x Ur sh ~ + Ar sh (|^ - (22) 
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In all the expressions deduced 



1 



n = 



R 0 cos B 0 



Now we proceed to derive the magnitude of the error 
arising from the fact that our solution is somewhat approx- 
imate as compared to the solution of the complete problem 
concerning the motion of the surface. 

Let us turn to figure 12, which represents a diagram 
of the nose of the dirigible straight model, the contour of 
which is given hy the equation y x = f(x x ). 

Let us take on the contour any point A, the trajec- 
tory of which during the circular motion of the dirigible 
coincides at a given moment with the direction of its 
speed u. The angle "between the direction of speed and the 
tangent to the contour we designate "by 8. The angle be- 
tween the tangent and the cross section AD we designate 
by n. Then, according to formula (9) and figure 9, we 
find that 



Now, in accordance with the solution given above, let 
us curve our model and designate the angle formed by the 
tangent to its contour, plotted through point A, and the 
direction of velocity in the tunnel by 6 2 . 

Analogously to figure 12, let us designate the angle, 
formed by the tangent to the contour and the tangent to 
the curved cross section AD, by |JL i , shown in figure 
13, which represents a part of the contour of the curved 
model y = f(x), the curved axis z ' D = f x (x) and the 
section AD. Then, in accordance with the principle for 
constructing cross sections, we get: 



8 + |A = tt — d 



(23) 



Pi + Hi = TT - ft 



(24) 



comparing formulas (23) and (24), we get 



Let us find the value of (J. i and compare it with p,. 
With these values we shall be able to estimate the differ- 
ence between 8i and B. From figure 12, we see that 
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u, = — - arc tan — — (25) 
2 dxi 

On figure 13 let us plot the following: Through point 
A let us draw the tangent MA to the cross section^ AD. 
Let us draw AN perpendicular to AM. Let us now give 
the argument xi the small increment Axi = DD ! . Then 
section DA will assume the position D 'A 1 . From point 
A' let us drop the perpendicular A»F on the prolonga- 
tion of the direction of velocity v and the perpendicu- 
lar A'E on the straight line AN. Let us designate by 
G the point of intersection of A'E and AN. Let us 
designate the angle "between the prolongation of the chord 
A 'A and the direction of velocity "by Pi 1 , and the an- 
gle between the same prolongation and the tangent MA "by 
Hi 1 . We readily see that 

p^t = Tr - $ - p x ' = arc tan j~ - ? 1 1 



But 



A»F 



Pi 1 = Z_ A'AF = arc tan ^ 



Noting that AF is an increment of x, and A'E is an 
increment of y for the curve of the curved contour y = 
f(x), we may write the preceding equation as follows: 



Then 



Ay 

0, ' = arc tan — — 
1 Ax 



111* = arc tan - arc tan ^ (26). 



dx 

hut L- A 'AF = l__ A ' AE + /_ GAF. 

Hence arc tan ^ = arc tan + /_ GAF 



"but since 



l_ GAF = tt - d = arc tan 4~ - ? 

dx e 



then the preceding equation may he rewritten as follows: 

Av i A'E . i. dz' tt 

arc tan = arc tan -r=— + arc tan -r — - -5 
Ax AE dx 2 
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Ay- 
Substituting the obtained value arc tan ^ in ex- 
pression (26), we find: 

, tt . A'E 

H X I = g - arc tan — 

It will be readily seen that, when Axj. infinitely tends 
toward zero , 



Ax 



lim = lim fn - arc tan fi^) 

i — * 0 ^ 1 Axj.— 0 V2 AE / 



Let us designate AE = Axi ' and rewrite the last equa- 
tion as follows: 

. . [tt /Ay i A'E Ax x \"| 

^ " Ax x ii. m o [2 - arc tan tsr Axi«;j 

where Ay x is the increment of the arc AD; that is, to 
say , 

Ay x = A'D - AD 

It is readily observed that when Ax x — * 0 , the ratio 

A'E Ax x 

g— tends toward unity, the expression "£5^* tends to- 
ward the value of the derivative \; in this way, 

dx ! 1 

■H - I - arc t« ^) (27) 

By the aid of the obtained expression (27) and the value 
|j, from formula (25), we find that the relation of the 
tangents Ho. and p, will be: 



tan 




CO t 


arc tan ( - — 


dXj. \ 






0. 


tan 






cot arc tan 


dx x J 





dx x ' 

dx7~" 



(28) 



The derivative — — obtained characterizes that de- 

dx x 

gree of contraction or elongation which is experienced by 
the lee side and the windward side of the model during the 
process of curving. 
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We can prove that the curved axis of the model. ap- 
proximates the arc of the circle with the radius p = 
R Q cos (3 0 , and the cross sections are almost plane. Mak- 
ing use of these circumstances for ascertaining the order 
of the numerical value of the ratio Hi/|J.> we find that 
for any point on the axis 

d *i 1 ± y\ 

dx x ™ R 0 cos 3o 



where y is the local radius of the cross section. 

It will "be readily seen that for the stern point and 
the nose point, where y x = 0, we will have y, = and, 

dxi ' 

consequently, 3 = 3 X . ^— — is farthest from unity in the 

vicinity of amidships. Let us find this value, assuming 
that R 0 cos 3 0 ~ 2.5 L and y « -| L (which corresponds 
to the mean given data for the model) J we obtain — » 

0.96; that is to say, the difference "between tan jj, and 
tan p, 1 amounts to ~ 4 percent. But, since in the vi- 
cinity of amidships "both \i and [x x approximate , the 

4 percent difference in the tangents creates a difference 
in the angles amounting to a few tenths of one percent. 
In the parts of the zones located "between the nose and 
amidships and "between amidships and the stern, the differ- 
ence "between 3 and 3 X is greater (due to the fact that 
in these zones the angles (i and Pj differ from tt/2 
more than in the case of amidships). However, according 
to our calculations, the difference "between 3 and 3 X 

nowhere exceeds from 2 to 2.5 percent. In this way, the 
insignificance of the error resulting from the abovemen- 
tioned alteration in the lengths of the arcs of the merid- 
ional contour has "been demonstrated. 

3. METHODS OP DIRECT DESIGNING OP CURVED MODEL FORMS 



There are two possible methods: (1) A very exact 
method, which takes into account the curvature of the ■ 
cross sections, and (2) a simplified method, in which 
the cross sections of the model are assumed to he plane 
and normal to the catenary of the axi s. 
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The axis of the curved model is constructed in the 
same way in each method. 

(1) Exact Method 

a ) Constr u ction of the curve d a xis ." To "begin with, 
it is necessary to find the radius of turn with which it 
is desired to carry out the test. In addition, we must 
know the angle of attack 8 0 formed during the rotation 
with the given radius. 

Knowing R Q , B Q may "be obtained either according to 
the previously constructed nomogram of the straight line 
~ = f(8 , 6) for the given dirigible, or from free- 

flight tests, or "by the aid of the statistical method. 
The last method is the simplest, because for the straight 

_ ~ R S 

line -~ = f(B 0 ), it gives the angle coefficient -Q—SSL - 

k, where a is the arm of the tail surface center of 
pressure. According to a whole series of foreign inves- 
tigations, this coefficient preserves sufficiently well 
its magnitude in the cs.se of quite a large number of air- 
ships. With k given, the angle of turn is calculated 
according to the formula 

Po = k (29) 

Here 8 0 is obtained in radians.** 

Knowing R 0 and 8 0 , it is necessary to find the 
distance from the nose of the model to point B, which is 
taken as the origin of the coordinates of the catenary. 

Designating by X x the distance from the nose of the 

model to the center of volumetric displacement, wo find, 
according to figure 9, that the distanco sought is: 

x iB = - x lQ = Xi - R 0 sin Po ( 3 °) 

♦See our paper "Determination of Radii and Angles of Turn 
of the Airship at Various Rudder Angles of Devia- 
tion." 

**In the experimental part of this paper we shall see that 
R Q and 8 Q can be assumed entirely arbitrarily. 
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If we obtain x^ in the negative, it means that in figure 

10 the nose point will lie to the right of the origin of 
the coordinates. 

For constructing the catenary of the axis z' = f(x), 
it is necessary, for calculating z* according to formu- 
la (7), to find the magnitudes x. The simplest way to 
find the magnitudes x, is to determine them according 
to formula (8), where x x is the distance along the arc 
of the catenary from point B. 

In this manner we can obtain on the curved axis of 
the model a series of points, the distance of which from 
point B and, consequently from the nose of the model 
along the arc is known. This is very convenient, since, 
in the profile of the noncurved model , the ordinates y x 
of the contour of the meridional section, as well as the 
position of the power cars, the tail surfaces, and other 
details are determined "by the distance from the nose. 

b) Now we have to proceed to the construction of the 
lines of the cross sections. These lines are construct- 
ed for each desired soction according to formula (16), 
where the magnitudes Xp_ and x 1 for each point of the 

curved axis are known from the previous construction. The 
moving coordinate x we may choose according to our own 
judgment. Then we have to find on the constructed curves 
points A and B (fig. 11), "belonging to the curve of the 
meridional contour; for these points the ordinates y x 

are known from the straight model. These constructions 
are made according to the formulas (21) and (22), which 
give the abscissas of the points A and B from the origin 
of the coordinates. Figure 14, showing this construction, 
gives an idea of how the net of the . cat enari es of the axis 
and of the cross sections looks. 

(2) Simplified Method 

The simplification concerns only the construction of 
the cross sections, which are assumed to "be plane and nor- 
mal to the curved axis. That such an assumption can he 
made, i.e., to substitute for segments of the catenaries 
AB a straight line, is shown by the following consider- 
ation. 

Let us find the radii of curvature of the catenaries 
of the cross sections at their intersection points with the 
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axis (that is, to say, when x = Xjj), The radius of curva- 
ture for the flat-topped curve is expressed: 

,3/2 



d 2 z» 
dx 2 

Taking the expression dz'/dx according to formula 
(18) when x = xp, we find: 



dzi = 1_ 

dx nx x 

Differentiating formula (17), we have: 

d 2 z ' 1 /x - x-r, 1 \ 

-—5- = — ch ( £ ~ Ar sh } 

dx" x x \ xi nx a / 

which, when x = x^ , gives: 



d z' = 1_ L + TT^T 2 " 
dx s Xi V W x > 

substituting the expressions obtained in the formula for 
p , we find that : 

3/2 



1 

= Xj + --r- 



Substituting n = --„ we have: 

R 0 cos 6 0 ' 

R 0 3 cos 2 6 0 
p = Xl + 2 

x 1 

let us show that in all the cases the value of p is 
very large. From the last expression, we can obtain: 

p ■ 0 + @o* £) 1 — >°] 1 

where L is the length of the model. 
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When x x = 0 p = «; that is to say, the cross sec- 
tion is a plane. When Xj increases, the magnitude p 

decreases. p will have the smallest value with the maxi- 
mum possible « L. On page 21 of this report it was 

shown that the minimum radius is equal to from 2 to 3 L - 
on the average, 2.5 L. The angle of attack, corresponding 
to the given E 0 , is usually about 7 to 9°, and conse- 
quently, on 'the average, cos 3 0 = cos 8° = 0.99 and 
cos 3 B 0 = 0.978. 

Substituting the values indicated in the expression 
found, we get the smallest possible radius of curvature of 
the cross-section diameter: 

p . <* 71 
r mm 

Such a radius is obtained at the tail of the model, 
where the diameter of the cross section is very small, for 
which reason it is entirely possible to regard the cross 
section as a plane. 

In the midship section, that is to say, when y x = max 
and x x « |j. we get : 

p a 21 L 

r mean 

This shows us that it is entirely possible to make 
the cro ss . sections plane. Besides, they preserve their 
circular form, which is a great convenience in preparing 
the model. 

Now let us proceed to the method of constructing the 
plane sections. It was shown above that these sections 
must be normal to the curved axis. Let us find the equa- 
tion of the normal at any point on the curved axis. The 
equation of the normal is generally written as follows: 

St (t - «») -g oi- - ° 

Here P is a function in an implicit form, | and 
T) are the moving coordinates of the normal line sought , 
and x-p and z 'j, are the coordinat es ' of the point on the 
curve, through which the normal passes. 
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According to equation (7) also F = z'j) - ^ ch nx D + 
= 0, whence 

= 1 and SI - - S h ni B 

dz' 8x JJ 

Substituting in the equation of the normal the values of the 

derivatives and taking, according to formula (8), sh nxj) = 

nx, . we shall have: 
1 D 

| - x D + nx ir (Tl - z« D ) = 0 

hence any abscissa of the normal will "be equal to: 

x x (T) - z D ) 

i = x D - nx lD (Tl - z» D ) = x D - -^"-~-p~ (3D 

It is most convenient to construct as follows (fig. 
15): In the diagram there is plotted on the curved axis 
point D with the coordinates x D , z'j), through which 
point it is necessary to draw the normal. If we take for 
the second point on the normal its intersection with the 
axis x (that is to say, if we assume in formula (31), 
Tj = 0 ) , this method is not very precise when the normals 
are constructed near the nose of the model, where the line 
of the axis almost coincides with the axis x. Therefore, 
it is more advantageous to find the intersection of the 
normal with some straight line parallel to the axis x 
and located at a sufficiently large arbitrary distance 
BM = z from i t . Then, when T) = Z , equation (31) gives 
the distance £z, from the axis z', of the points N 



x 1D (Z - zt D ) 

Pz = x-n - — -r [61a. ) 

& u R 0 cos 3 0 



In order to find the points A and B, belonging to 
the contour, we simply have to lay off on the lines of the 
normals from the points D the magnitudes j x . 

It is self-evident that also the tail surfaces of the 
airship must be somewhat deformed. Both vertical stabil- 
izers must be bent along the segment of the catenary of 
the axis at the place where they are installed. But the 
horizontal stabilizers are constructed in the coordinate 
net formed (in the simplified case) by means of catenaries; 
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equidistant from the line-' of the axis, and by means of the 
normals to the catenaries. In this way, the right stabil- 
izer comes out shorter than the normal, and the left sta- 
bilizer i& longer than the normal. But the. dimensions 
laid off on the normals remain unchanged. 

The model is installed in the tunnel in such a way 
that the center of volume of the model lies on the axis of 
the tunnel, and the axis x (fig. 10) is placed parallel 
to the axis of the tunnel at a distance z, c > which is 

easily determined according to formulas (7) and (8), if 
we know the position of the center of volume on the axis 
of the model. 



4. THE VELOCITY GRADIENT ACROSS THE TUHNEL 



It was shown above that, in order to maintain the 
similarity of circular flight, it is necessary to create 
in the tunnel, when testing the curved model, a velocity 
varying, across the tunnel, according to a linear law. 
That this law is linear is evident from figure 9. At the 
moment represented on the diagram, the airship turns at a 
constant angle of speed u>y. Then the linear speed of its 
center of volume (point C) will be: 



V, 



and the linear speed of any point on the axis, for exam~ 
pie, point D, will be: 



V = 



Wy * 



Subtracting one equation from the other, we get: 

Wy (R - R 0 ) 



hence 



\3 0 + UUy (R - R 0 ) 



When passing from curvilinear flow to rectilinear flow, 
it is evident from figure 10 that for point D 



R ~ R 0 = z> - z« c = z 



where z denotes the variable ordinate, read off from the 
axis of the tunnel crosswise to the right and to the left. 
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In this manner 

u = v 0 + UJy z (32) 

The expression obtained, gives the distribution of 
velocity along the section of the tunnel. As is evident, 
this law is linear. Substituting for ujy its expression 
in terms of the radius of rotation and the circumferen- 
tial velocity, in the given case, the velocity of the cir- 

cular flight of the airshit>; that is to say, cuy = — , we 



shall have: 



V - Vo + g . = v 0 (l + £) 



(33) 



Hence, the constant velocity gradient along the section of 
the tunnel 



9u _ = 
3Z " R 0 



"7 



that is to say, numerically it is equal to the angular ve- 
locity. 

The di stribution of velocity obtained has a flow which 
is vortical at every point. It is interesting to note 
that the velocity circulation along any closed contour in 
such a flow is equal to r = S ' u>y , where S l is the area 
confined by the contour. In this manner, the intensity of 
the vortices of such a flow is obtained as X = u>y , that 
is, half as large as in the case of circular flow. Conse- 
quently, (see p. 4), in testing the curved model, we are, 
in form, even closer to full-scale tests than in the case 
of testing in circular flow. 

The law, expressed by equation (33), is diagrammatic- 
ally represented in figure 10. In order to have an idea 
of the magnitude of the velocity gradient in that portion 
of the cross section of the tunnel where the curved model 
is installed for testing, we shall calculate the relation 
of the velocity at the boundary of this portion to the ve- 
locity along the axis of the tunnel. From equation (33), 
we get: 
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Analogously to the preceding, let us assume: 

R 0 * 2.5 L 

and the entire width, of the zone, in which it is required 
to create the gradient, 

2z = 0,5 L 

whence z = 0.25 L 



Such a width of the zone is entirely sufficient in or- 
der to provide ample room for any curved model. Hence, 

JL « i + gx|5 L „ lml 
V 0 2.5 L 

i.e., the difference between the velocity at the "boundary 
of the portion and the velocity along the axis of the tun- 
nel amounts to about 10 percent. The methods of creating 
the velocity gradient and the method accepted by us will 
be discussed in the experimental part of this paper. 

5. METHOD OF DETERMINING- ROTARY DERIVATIVES PROM 
WIND-TUNNEL TESTS OE THE CURVED MODEL 



We have already mentioned (p. 3) that, in directing 
the air flow against the curved model in the tunnel, we 
should obtain the forces and moments acting on the model 
during circular motion. Let us denote the force of the 
drag during circular flight by ZX, the. lateral aerody- 
namic force during circular flight by ZZ lf and the mo- 
ment of yawing (around the axis y) during circular flight 
by EM. 

Then in the very general case, any of these quanti- 
ties for any one airship will be a function of the angu- 
lar velocity, of the linear velocity, of the angle of at- 
tack P 0 and of the angle of rudder deviation; i.e., for 
example, the force 

SZx = f(u) y , v 0 , 0 O , S) 

Assuming that the motion is effected with v 0 = const., 
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P 0 = const., and 8 = const., "but with variable Wyt* 
we may write: 

SZj = Cp(U)y) 

Resolving the function obtained according to the 
Mac Laurin system which, as is well known, is as follows: 



SZi = cp(0) + ^ qp'(o) + cp"(0) + ^ <p'" (0) + ... 

we find 
ZZ^ = (ZZ 1 ) 



Wy' 



ivy- 



ra(Ez x ) 



Oly 

+ — 

(JUy =0 1J L 3'Uly 



Wy 



IDy 



= 0 



2J 



■a^EZi) 



(Uy 



| + . . . 

-Jojy=0 



T,Z 1 is transformed into the 



But when (i)y = 0 the force 

force Z x , acting on the model in rectilinear motion. 
Thus , 



az. 



^ ujy a rf z x 

«Zi = Zi + ^y + — SSSp + 



= z x + z. 



uu 



(34) 



where "by Z 
series, containing ujy 



x we designated the sum o f all the terms of the 



As will be seen, the full force, acting in circular 
flight, is the algebraic sum of the effects of the linear 
and rotary displacements. 

Analogously, resolving in the series EX and EM, 
we finally find: 



EX 


- X + Xyj 








EZ X 


= z, + z. 






(35) 


EM 


= My + M^ 


J 






the 


terms X^, 


Z 


, , and M„,: 






= EX - X 


- 






Zi 


— £ 2 ^ Z j 






(36) 




= EM - My 


_ 







♦Naturally, such a motion appears forced (i.e., the forces 
and moments acting on the airship are not in equilibrium ) 
and corresponds either to motions occurring at different 
times along circles of different radii, or to the motion 
along the noncircular trajectory. 
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which, as will he seen, are defined as the difference "be- 
tween the results of tests with the curved model and the 
straight model. Naturally, in addition to this, "both for 
the curved and the straight models, we must have 6 0 = 
const., v Q = const., and 8 = const., which was the ini- 
tial assumption in deriving the expressions (35) and (36). 
The forces and the moment X^, Z 1{1J , and M w , similarly to 

the usual aerodynamic factors, can he put in the form of 
nondimensional coefficients. 

Thus, as will be seen (considering the indicated non- 
dimensional coefficients as functions of the angular ve- 
locity), the method of curved models permits us, for the 
solution of many dynamic problems, to dispense with the 
measurement of the rotary derivative. 

But, paying tribute to the generally accepted method 
of studying curvilinear flight by the aid of rotary deriva- 
tives, we shall find the latter. 

As is well known, by rotary derivatives we understand 
the increment of the rotary moment or force per unit of 
angular velocity, that is to say, 

Of course, it would have been more exact to write in 
the left members of the equations, not dX , 9Zi and SMy , 

but dEX, d£Zj, dEM. But the designations for the deriv- 
atives here given are simpler and have already been used 
by us in our earlier work. 

The rotary derivatives are generally referred to the 
unit of speed of flight. Writing in place of X^, Z^, and 

Myj their values from (36), we divide the right and left 
members of the preceding equations by Vq ! 
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Substituting in the right members u) y = ~, we get: 



I- 1 - -5-1 

[ 



'EZ X 



1 9M y - 
u 0 Bujy. 



" »o 2 ) 



In the parentheses of the right members were obtained 
the relations of the aerodynamic factors to the squares of 
the velocities, which it is customary to designate by the 
symbol R with the subscripts x, z x , and m y . Introduc- 
ing the designations 



H X - R 

u 2 " EX ! 



EZ i - R 



v 



EM 

2 = R 



EM" 



we finally get 



R o ( R ZX - R X) 
R o (HeZj - R Z 1 ) 



"_1_ dX_~| _ 

PL iiil = 



(38) 



The expressions obtained enable us to find the rotary 
derivatives from wind-tunnel tests of curved models by 
measuring on the balance and on the moment apparatus the 
aerodynamic force and moment when the model is installed 
as shown above. The factors R x , R Zl . and R my , in the 
case of the angle 0 O , are known from the normal tests. 

It is clear that, strictly speaking, there is no 
sense whatever in measuring the forces and moments for the 
curved model at its various angles of attack with respect 
to the flow (as is done with straight models), since for 
each angle 3 Q there is a corresponding, degree of curva- 
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ture of the model;* this may be seen, for example, from 
formula (?) . ... 

However, we must .remember that the value we obtain 
for the rotary derivative of the drag is a true value 
only for the model, since, being determined by the coef- fit 
ficients of. drag, the value will depend to a large ex- ^ 
tent on the Reynolds Number; this cannot be said of the 
rotary derivatives of the lateral force and moment. 

In this manner, the method of curved models gives us 
(as will be seen from the expression (37)) the values of 
the rotary derivatives - in the first place, independently 
one from the other and, in the second place, taking into 
account all the terms of the series (34); that is to say, 
the generally accepted assumption that the rotary deriva- 
tives are independent of the angular velocity not only was 
nowhere assumed, but may be verified experimentally. 

This constitutes the colossal advantage of the method 
analyzed, ~- — 

6. VARIANTS OP THE APPLICATION OP THE METHOD 

OF CURVED MODELS 



Now, let us see in what variants can the method of 
curved models be applied in the aerodynamical investiga- 
tions of an airship. 

We have already shown above that the values R 0 and 
f3 0 , which are necessary for constructing the curved model 
profile, can be taken from the nomogram of the straight 
line : 

- 1 - = f(P 0 ,S) 

constructed according to the method set down in our work 
on radii of turn. In such a case the testing of the curved 
model will be an excellent verification for the values of 
the rotary derivatives, determined by the method of small 
oscillations* 



*In the experimental part of this paper it will be shown 
that this strict condition need not be adhered to, at the 
expense of a rather small error, thus extending the limits 
of application of the method of curved models. 
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Thus, the first variant in the application of the 
curved model is the verification of the values which we 
have for the rotary derivatives. Let us examine this va- 
riant in detail. 

The initial data (R Q and j3 Q ) for the construction 

of the curved model profile are taken from graph — = 

E 0 

f ( 9 0 t 6 ) • This graph, constructed on the oasis of the ex- 
periments with the model of the 7-2 airship, is repre- 
sented in figure 16. It shows the intersection of two 
families of curves, one of which gives the value of l/R-o 
in the presence of the condition of equilibrium of the 
forces only, acting on the airship during rotation.* The 
other family gives the values of l/^o ^ n the presence 
of the condition of equilibrium of the moments only. 

Naturally, the line (straight line), joining the 
points of intersection of the corresponding curves of the 
first and second families, gives the values of ~- = 

K 0 

f (3 0 »$) when there is equilibrium both of the forces 
and of the moments; i.e., it gives the values of R 0 , 9 0 , 
and 5, corresponding to the flight along circles which 
is established. 

Now, let us assume that we are preparing six curved 
models, each of which will be constructed by the aid of 
E 0 and & D obtained from the points of intersection 

of the corresponding curves of the first and second fami- 
lies. During the tests in the tunnel, the rudder of each 
of these models should be deviated at that ' angle 6, to 
which, according to figure 16, correspond the points of 
intersection of the curves of both families. 

In the result of wind-tunnel tests we obtain the 
force EZj = Zi + Zi^, which will be equal numerically 

to the centrifugal force (from the condition of equilib- 
rium). According to the force ZZ^, we find the rotary 
derivative (formula (37)). 

Thus, the aerodynamic moment of each of the curved 
models, when the models are properly installed in the tun- 
nel, will be equal to zero since, under the conditions of 
equilibrium of moments, the aerodynamic moment is extin- 



*The centrifugal force also participates in the equilib- 
rium. 
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guished "by the rotary effect, and according to the third 
formula of the system (38) the rotary derivative of the 



In this manner, the very result . R^M = 0 shows the 
accuracy of the rotary derivative of the moment, deter- 
mined by the method of damped oscillations. In order to 
obtain therefrom the generally accepted relations of the 
rotary derivatives to the angle of attack, it is necessary 
to construct the rotary derivatives, obtained from the 
wind-tunnel tests of the entire series of curved models, 
as a function of the angle 8 0 (which, for each model, 
was a basic factor in the construction of its profile). 

If we extrapolate the obtained graphs of the rotary 
derivatives up to 8 0 = 0 , we obtain the values of the 
rotary derivatives which are necessary for calculating the 
criteria of dynamic stability. 

Now let us see how the method of curved models can be 
applied as an independent method. 

Let us assume that for a given whole series of values 
of R 0 and B 0 , we construct as many models as the com- 
binations that can be made from R 0 and 8 0 . Naturally, 
each of the models will have its own degree of curvature. 
Each model must be tested on EM and EZ X for different 

values of the rudder angle of deviation 6. After this, 
grouping the obtained values EM and EZi for the vari- 
ous 6 and plotting them as a function of R 0 for each 
angle P Q , it is necessary, by means of graphical interpo- 
lation, to select on the curves EZj. = f (R 0 ,8) those val- 
ues of R Q and 0 O for which, from the, condition of equi- 
librium of forcos, EZj, is equal to the .centrifugal force 

-t£°— cos 6 0 (by intersecting the net of curves EZ X = 

m V s 

t (R 0 ,6 0 ) with the net of hyperbolas _ r~ cos Po = qp(Ro))» 

Thus, having those values of ' R 0 and B 0 with which, 
when 6 is given, equilibrium of the forces is. obtained, 
we obtain that family of curves which in figure 16 charac- 
terizes this equilibrium. 
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Prom the curves ZM (R 0 ,8) we select those values of 
R 0 and 8 0 for which, according to the condition of equi- 
librium of moments, EM = 0, 

In this manner we obtain the second family of curves 
(fig. 16), as a result of the intersection of which with 

the first family we obtain the sought nomogram — = 

K o 

f(0o» 8 )» which is determined with incomparably greater 

accuracy (that is to say, taking into account the influ- 
ence of all the variables), than with the method set forth 
in our work on radii of turn. 



It is obvious that we can also obtain the graphs of 
the change of the rotary effects of the moment, of the 
lateral force and of the drag as a function of the angle 
of attack, of the rudder angle of deviation and of the an- 
gular velocity. 



Besides, by testing only curved hulls of airships 
(without tail surfaces) and working up the data by the 
aid of methods analogous to the one described, we can ob- 
tain., simultaneously with the very interesting relation, 



i a it 



y 

v 0 a~ai v 



= f(Po) 



Trull 



(which it is extremely difficult to obtain in the tunnel 
with the method of small oscillations) as well as the val- 
ues : 



[ 



_1_ BZ^ 

u 0 3UUy 



= fO 0 > 
hull 



which, as far as we know, have not been determined any- 
where experimentally with sufficient accuracy.* 



*The English, in their work, calculate for all airship 
hulls in general, 



azi 



r- - 

LUo -au) y j hull 



0.1 m 



where m is the mass of the airship; this requires exper- 
mental -or oof. 
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It is self-evident that such a complete investigation 
of curved models of airships is profitable only if we have 
one model which can he curved to the desired profile. The 
question of the possibility of producing- such models (elas- 
tic) will he taken up next year by the Air Section of the 
Experimental Aerodynamic Department of the Central Aeror* 
hydrodynamic Institute. 

By arranging the pressure tubes of the curved models, 
we can find the distribution of pressure in flight along 
circles. 

By the way, it is interesting to note that, by the 
aid of curved model tests on the distribution of pressure, 
we can find experimentally the influence on the final re- 
sult of that altered length of the arc of the meridional 
contour which is obtained when the curvature is effected* 
According to our opinion, it is necessary for this purpose 
to integrate the distribution of pressure along the sur- 
face of the curved model and compare the obtained result 
with the integral of the same distribution of pressure, 
but with the corresponding points located along the sur- 
face of the straight model. 

Thus, we think that the method of curved models here- 
in expounded, when used on the very widest scale, can and 
must become the most exhaustive method of studying the dy- 
namics of curvilinear flight and the dynamic stability of 
airships. 



Translation by Translation Section, 
Office of Naval Intelligence, 
Navy Department, 
Blum a Karp, 
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Figure 1. 
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Figure 2. 
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Figure 4. 
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Figure 5. 
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Figure 7. 





Figure 15. 
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